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The filamentation instability occurring when a non relativistic electron beam passes through a 
quantum magnetized plasma is investigated by means of a cold quantum magnetohydrodynamic 

QO ' model. It is proved that the instability can be completely suppressed by quantum effects if and 

only if a finite magnetic field is present. A dimensionless parameter is identified which measures 
the strength of quantum effects. Strong quantum effects allow for a much smaller magnetic field to 

£N| ■ suppress the instability than in the classical regime. 
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I. INTRODUCTION 



The development of quantum hydrodynamic and magnetohydrodynamic equations [l|, Q made it possible to quickly 
evaluate quantum effects connected to the physics of microelectronic devices and laser plasmas interaction (see Ref. 
Q and references therein) . Plasma physics has also gained from these progresses as quantum effects appear in Fusion 
settings or Astrophysics. The behavior of waves in quantum plasmas [J 0, H, B 0] magnetized or not, as well as 
turbulence is such environments Q has thus received attention. Another very classical topic of plasma physics, 
namely plasma instabilities, needs to be revisited from the quantum point of view. The quantum theory of the 
two-stream instability has already been developed (9l. ITo| while quantum effects on the filamentation instability were 
recently evaluated for a non magnetized plasma. Due to the importance of magnetized plasmas, especially in 



astrophysics, we devote the present paper to the evaluation of quantum effects on the filamentation instability in such 
setting. Since the relativistic quantum magnetohydrodynamic equations are yet to be defined, the present analysis 
restricts to the non-relativistic regime. On the other hand, we do not make any approximation on the beam density 
so that present theory remains valid even when the beam density equals the plasma electronic one. 
Oh The paper is structured as follow: we start explaining the formalism and derive the dispersion equation. We then 
turn to the investigation of the marginal stability and derive some exact relations satisfied in this case. We finally 
I , study the maximum growth rate and the most unstable wave vector before we reach our conclusions. 
■ Let us then consider an infinite and homogenous cold non-relativistic electron beam of velocity Vf,z and density n& 
entering a cold plasma along the guiding magnetic field Bo = Bqz. The plasma has the electronic density n p and 
ions form a fixed neutralizing background of density n\, + n p . The beam prompts a return current in the plasma with 
velocity V p z such as n p V p = n^Vb- We use the fluid conservation equations for the beam (j = b) and the plasma 
U =P), 



^ + V-(n i v i ) = (1) 

o ■ 

^ ' and the force equation in the presence of the static magnetic field Bo with a Bohm potential term [2j , 

+ Vj • V)v, = —3 ( E + + — rV —*=l , (2) 



<9v, — o v, x B 

53 ■ -^ + (v ^ V ^ = -^( E+ ^-r 2m 2 

where q > and m are the charge and mass of the electron, rij the density of species j, pj its momentum, and B 
equals Bo plus the induced magnetic field. We now study the response of the system to density perturbations with 
k _L Vfc, varying like exp(ikr — tot) with k = fcx, and linearize the equations above. With the subscripts and 1 
denoting the equilibrium and perturbed quantities respectively, the linearized conservation equation (TT]) yields 

k ■ v 7 "i 

nji=rijo : , (3) 

u> - k • Vj-o 



'Electronic address: antoincclaudc.bret@uclm.es 



2 



and the force equation ([2]) gives 



q ( v j0 x Bi + Vji x B \ .hk 2 riji 



m \ c / 4m 2 rijo 

From the linearized equations above, we derive the perturbed density and velocity fields in terms of Ei and Bi and 
eventually express the current through, 

J = q ^2 rijoVji + rijiVjo. (4) 

j=P,b 

Finally, we express Bi in terms of Ei through B x = (e/aj)k x Ei and close the system inserting the current expression 
in a combination of Maxwell Ampere and Faraday's equations, 

C 2 4?7T 

— kx (kxEi) + Ei + J = 0^T(Ei) = 0. (5) 

The tensor T has here been calculated symbolically using an adapted version of the Mathematica Notebook described 
in Ref. [la . It takes the form 

/Tn T? , \ 

(6) 
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where the superscript * refers to the complex conjugate and 
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T 12 - 1 ( x 2 _ p2 \ S^4' ( 7 ) 

in terms of 



x = — , Z = , p= — , a = — , \l B = , (8) 

Up u>p c n p mcujp 

where uj p is the electronic plasmas frequency. Quantum effects appear to be measured through a parameter previously 
highlighted [II [13, ED], 



6 = f> 



(9) 



Numerically, 

6 C = 1.3 x 10" 33 n p [cm" 3 ], (10) 

so that this parameter will hardly be larger than 1, even when dealing with the densest space plasmas. 

Let us finally emphasized a point regarding the beam to plasma density ratio a defined in Eq. (|8[). If the ground 
state which stability is investigated consisted of the plasma only, the beam representing the perturbation, then this 
parameter would have to remain much smaller than 1 within the framework of a linear response theory. In turns out 
that the dispersion equation which has just been derived is the dispersion equation of the beam+plasma system. The 
perturbed ground state is therefore the sum of the beam and the plasma. It is thus perfectly possible to investigate 
the linear response of the whole system even when a — 1 so that we do need to make any assumption regarding this 
parameter. 
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FIG. 1: Classical (non-quantum) growth rate of the filamentation instability in terms of the reduced wave vector Z without 
(a) and with (b) magnetic field. Curve (c) includes quantum effects. Parameters are a = 0.1, /3 = 0.1 for (a,b,c), Q.b = 0.03 
for (a,b) and Q c = 1.3 x 10 -7 for (c), which corresponds to the plasma density n p — 10 26 cm -3 . 



II. CLASSICAL MAGNETIZED PLASMA 



Before we turn to the quantum case, let us quickly remind some basic features of the cold magnetized filamentation 
instabilities [l3[ in the classical (non-quantum) regime. To this extent, the dispersion equation, which is just the 
determinant of the tensor we just defined, is solved numerically and Figure [1] displays the growth rates obtained with 
and without the magnetic field (curves a and b). The stabilizing effect of the magnetic field is twofold. On one hand, 
the smallest unstable wave vector switches from Z — to 

Zl = - (11) 
y/a(l + a)[3 2 - Q 2 B 

On the other hand, the growth rate saturation value doo for large Z is lower with 

^ = yJa(l + a)P*-n 2 B , (12) 

which vanishes exactly for 



n B = n Bc ee /V«(i + a) (13) 

Noteworthily, this value of the magnetic field also makes the quantity Z\ diverge. The physical interpretation of this 
threshold is simple as /3y a(l + a) is just the maximum growth rate of the instability in the non-magnetized case 
[p4| . Filamentation instability is thus inhibited when the electron response to the magnetic field is quicker. 



III. QUANTUM MAGNETIZED PLASMA 
A. Marginal stability analysis 

Figure [TJ; displays the growth rate in terms of Z accounting for quantum effects. As in the non-magnetized case 
[ll| . quantum effects introduce a cut-off at large Z so that we now have to characteristic wave vectors Z\ and Z 2 
determining the instability range. Both of them can be investigated directly from the dispersion equation. Since the 
growth rate vanishes for these wave vector while the root yielding the filamentation instability has no real part, we 
can write 

det T(x = 0, Z = Z x ,-i) = 0. (14) 

It turns out that this equation can be simplified. After replacing Z 2 — > Z and eliminating x — as a double root of 
the dispersion equation, we find that the equation above is equivalent to P{Z)Q(Z) = with 



P{Z) = (e c Z 2 +p 4 n 2 B )(Z + (l+a)p 2 )-Za(l + a)p 6 
Q{Z) = Z/3 4 fi| + (Z + (l + a)/3 2 )(O c Z 2 + (l + a)/3 4 ). 



(15) 
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FIG. 2: Schematic representation of the functions F and G defined by Eqs. (I16|) . 



Every term of the second polynomial is clearly positive so that it yields only negative roots Z < 0, implying some 
complex wave vector. Because we seek real wave vectors, we can conclude that Z\ = {Z\) 2 and Z^ — {Z2) 2 are both 
zero's of P(Z). This function being a polynomial of the third order, it is possible to find the exact solutions. We 
nevertheless use some graphical method for a more intuitive approach. Let us then define 

F(Z) = (e c Z 2 + f3 4 n 2 B )(Z + (l + a){3 2 ), 

G(Z) = Za(l + a)p 6 , (16) 

so that P = is equivalent to F = G. F is a third order polynomial, monotonically increasing for Z > 0, and 
starting from F(0) = (1 + a)Q 2 B l3 e with an initial slope F'(0) = G is a first order monotonically increasing 

polynomial with G(0) = and slope a(l + a)(3 6 . We can now conduct the graphical analysis of the problem following 
the guidelines set by the schematic representation of F and G on Figure[2l When increasing Q,b or C , the curve G 
is not modified because neither Qb nor O c appear in its expression. Meanwhile, F(0) increases with Qb, and F(Z) 
increases all the more than O c and f2s are large. This allows us to draw the following conclusions: 

• In the absence of magnetic field, F(Q) = ^'(0) = while the previous analysis remains unchanged. The equation 
F = G thus has two positive solutions regardless of the other parameters. One solution is Z\ — 0, i.e. Z = 0, 
and we label the other Z2 > 0. We recover the existence of a quantum cut-off [llj at large wave vector, and 
prove here that the instability is never completely stabilized since Zi never vanishes. 

• For any finite magnetic field, one has L^O) > and L"(0) > 0, and the typical resulting situation is the one 
represented on Fig. [5] As long as F(0), or the growth of F(Z), are "not too high", the equation F = G has 
two positive roots Z1.2 = [Z\ t i) 2 . But it is obvious that as FIb or C increase, the too real roots become one 
before they vanish. We thus come to conclusion that the instability can be completely suppressed by quantum 
effects if, and only if, the system is magnetized, regardless of the strength of the magnetic field. It is graphically 
obvious that since an increase of both £Ib or 9 C contribute to the collapse of the two reals roots, the stabilization 
condition should result in a balance between these quantities. It should be possible to stabilize the system at 
low Qb with an high C , or vice versa. 

In the magnetized case, the instability is marginal when the equation F — G has one double root Z\ = Zi = Z^. 
Here, L stands for Last because Zl = \[Z~Z is eventually the last unstable wave vector before complete stabilization. 
If Zl is double root of F for the parameters defining the marginal stability, then for these very parameters F can be 
cast under the form F(Z) = (Z — a)(Z — Zl) 2 where a is the third root. By developing this last form and identifying 
the coefficients of the polynomial with the ones extracted from Eq. (|16p . we can write the following equations, 

a + 2Z L = -(l + a)p 2 , (17) 

2aZ L + Zl = + (18) 

aZl _ _(i±^k_ (19) 
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FIG. 3: (Color online) Values of O* and Q* B implicitly denned by Eqs. (|20I21[1 for various beam to plasma density ratios a 
and /3 — 0.1 (red bold curves) and /3 = 10 -3 (blue thin curves). Given a and f3, the system is marginally stable for parameters 
located on the corresponding curve, and stable above. The value of Qbc (see Eq. 1 13 p for /3 = 0.1 and a = 1 is 
represented by the horizontal dashed red line, and the oblique dashed curve corresponds to Eq. (J25J with the same a, f3. 



where the superscript * refers to the values at marginal stability. By eliminating a between the first and the second 
equation, one finds a second order equation for Zl which positive solution can be cast under the form, 



Then, eliminating a between the first and the third yields an implicit relation between Zl, 0* and £l B at marginal 
stability, 

(1 + a)/3 6 n*i = Z 2 L e* c (2Z L + (1 + a)(3 2 ). (21) 

Equations (|2QI2ip therefore define 0* and £l* B in terms of each other, and of the others parameters of the problem. 
The curves thus defined appear on Figure[3]for various a's and (3's. Parameters (0 C , fis) located above a given curve 
(<d*,£l* B ) define a completely stabilized system. 



IV. ANALYTICAL EXPRESSIONS FOR MARGINAL STABILITY 



A. Classical limit 



We observe on Fig. [3] that £l* B reaches a finite value when O* — > 0. This classical limit is obviously the marginal 
magnetic parameter Qbc given by Eq. (fT3|) . We thus assume a leading term in the development of £l* B for small 0* 
of the form Q* B = (1 — K<d*t)Q Bc . Inserting this expression in Eqs. (|20|21|) and expanding the results in series of 0*, 
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we find 



and 



i 1 - ^"^ : 173 ) "Be, (22) 



gV ^l + g) 2 / 3 ^/ 3 
2V30* 



In accordance with the classical case where the smallest unstable wave vector diverges for marginal stability (see Eqs. 
Ill|12p . the last unstable wave vector Z^ behaves like l/O* 1 ^ 6 in the weak quantum regime since = Z L . 



B. Strong quantum limit 

Having elucidated the weak quantum regime, we now turn to the strong quantum one. Figure [3] makes it clear that 
marginal stability behaves differently within each regime. In order to discuss this point, let us consider expression 
(l20f of Zi in terms of the marginal classical magnetic parameter VLb c - In the "large" 0* regime, the ratio under the 
square root becomes small compared to unity, and Fig. [3] shows that £l* B <C Qbc- Developing the square root, wc find 
directly 



20*(l + a)' 



(24) 



In this strongly quantum regime, the last unstable wave vector thus tends to zero like Inserting the former 

expression in Eq. (I21|) yields the magnetic parameter required to stabilize the system 

2(l + a)0*- 2A /0*' { ) 

This limit is plotted on Fig. |3Jfor f3 = 0.1 and a = 1 and perfectly fits the numerical evaluation for large 0*. 

It is now possible to exhibit the dimensionless parameter measuring the strength of quantum effects. The equation 
above indicates that fl B <C £Ibc if ^Bc "C 2(1 + cx)Q*, and the curves plotted on Fig. [3J demonstrate that a reduction 
of the stabilizing parameter is the signature of the strong quantum regime. Because we think here in terms of orders 
of magnitudes, we drop the 2(1 + a) factor and finally define 

A - , (26) 

as the parameter determining the strength of quantum effects. These are weak for A>1 and strong in the opposite 
limit A < 1. 



V. UNSTABLE SYSTEMS 



Having elucidated how the system can be completely stabilized by quantum magnetic effects, we now turn to 
unstable systems in order to investigate the growth rate of the instability and the most unstable wave vector for a 
given configuration. 



A. Maximum growth rate 



In the weakly quantum regime with A ^> 1, the maximum quantum growth rate is very close to its classical 
counterpart all the way down to complete stabilization which, as we just mentioned, occurs for similar magnetic 
parameters fis (see Eq. [22] above). Figure [4] present a plot of the maximum growth rates along the Z axis, in terms 
of fig in the classical and quantum cases for (3 — 0.1 and various a's. Parameters have been chosen to illustrate the 
present weak quantum regime. Such a system can thus be viewed as basically magnetized with some weak quantum 
effects, and stabilization mainly comes from the magnetic field. 
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FIG. 4: Maximum growth rates in terms of Ob in the classical and quantum cases for = 0.1 and a = 1CT 1 , 1(T 2 , 10" 3 . With 
Be = 10 -10 (n p = 7.6 x 10 22 cm -3 ), the parameter A given by Eq. (|26[l is always larger than 3 x 10 7 and the classical and 
quantum curves are hardly distinguishable although the quantum growth rate is a little bit smaller. 



When A <C 1 (strong quantum regime), stabilization is reached earlier with respect to Qb (see Fig. [3] and Eq. |25|) . 
Here, stabilization comes from a combination of quantum and magnetic effects, as indicated by the oblique slope of 
the curves in Fig. [3l We plot on Figure [5] the maximum growth rate in terms of Qb ■ We recognize the kind of curve 



obtained for a classical magnetized plasma S. 



2 

max 



B 



fl B with a "cut-off" magnetic parameter a/3 /2yQ* 



This is why we plotted together the numerical evaluation of the maximum growth rate together with the function, 



\ 



a(3 2 



- n 2 



(27) 



It can be checked that this function fits the result all the more than A is small. With Eq. ([23]) . we then come to 
the conclusion that as far as the maximum growth rate is concerned, strong quantum effects are equivalent to the 
substitution, 



n Bc & A 



Be 



2(1 + a)' 



(28) 



Because this new quantum cut-off is much smaller than the classical one, the maximum growth rate is reduced 
accordingly. 



B. Most unstable wave vector 



The most unstable wave vector is, together with the maximum growth rate, the most relevant information about 
the unstable system. In the classical case, the growth rate just saturates at large Z yielding a continuum of most 
unstable modes. But quantum effects stabilize the large Z modes, so that there is always one mode growing faster 
that the others. 

For systems near marginal stability, the last unstable wave vector Zl = y/ZE, given exactly by Eq. (|20|) . and in the 
weak and strong quantum limits by Eqs. (|23I24[) respectively, is by definition a very good approximation of this most 
unstable wave vector when replacing the marginal parameter 0* by its actual value O c . Indeed, we found numerically 
that expressions (I23I24[) are still quite accurate, even for systems far from stabilization. This can be understood from 
Fig. [5] on one hand, the most unstable wave vector for a given configuration is necessarily between Z\ and On 
the other hand, the last unstable wave vector belongs to the same interval because Z\ increases while Zi decreases 
as the system moves towards stabilization. For a typical situation such as the one represented on Fig. [TJ:, Z\ and 
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FIG. 5: Maximum growth rates in terms of f2s for /3 = 0.1 and C = 10 . The thin curves have been computed numerically, 
and the bold ones (when distinguishable from the thin) correspond to Eq. (|27|l . Agreement improves with A decreasing. 



Zi are eventually quite close to each other so that Zl, which is in between, cannot be far from the most unstable 
wave vector. With the parameters chosen for this plot, we find A = 2.5 x 10 5 indicating a weak quantum regime. We 
therefore turn to Eq. (|2"5|) and find the most unstable wave vector for Z ~ 0.4, which fits accurately what is observed. 



VI. DISCUSSION 



Quantum effects have been assessed with respect to the filamentation instability in a magnetized plasma. As far 
as the unstable wave vector range is concerned, magnetic effects set it a finite lower bound, while quantum effects 
introduce a cut-off at large A:. As a result, the unstable domain takes the form [ki, ^2] and can eventually vanish for 
some parameters configurations which were elucidated. 

We also found that the dimensionless parameter A = Qbc/Qc determines the strength of quantum effects. When 
A ^> 1, the instability can be described in classical terms, and eventually vanishes when increasing the magnetic field, 
while the unstable wave vector range shifts towards infinity. When quantum effects are strong, namely Ac 1, the 
instability still vanishes with the magnetic field, but the unstable wave vector range tends to zero. Furthermore, the 
magnetic field required to stabilize the system is divided by A/2(l + a)<l with respect to its classical value, so that 
filamentation can be suppressed by a much smaller magnetic field than in the non-quantum case. These results may 
have important consequences when dealing with dense space plasmas. 

Finally, it will be necessary to assess both relativistic effects, which tend to enlarge the instability domain while 
reducing the maximum growth rate [l5|, and kinetic effects which usually have a stabilizing effect (l6j . To this extent, 
relativistic quantum kinetic theory will be required, or the relativistic form of the quantum Euler equation ([2]) will 
have to be elaborated. As long as the theory implemented is non-relativistic, the magnetic stabilization level unraveled 
here should remain an upper stabilization bounds when kinetic effects are accounted for. In the classical relativistic 
regime, it has been demonstrated that the stabilizing magnetic field behaves like ^% where % is the relativistic 
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factor of the beam. Because this increase of the magnetic threshold eventually stems from the relativistic increase of 
the mass of the electrons, we can conjecture that the same factor ^/t^ will be found in the relativistic counterpart of 
Eq. (|25|) . but this shall need confirmation. 
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